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Abstract 

"(N : 

Two varieties of the universal stellar initial mass function (IMF) viz., the Kroupa and the Chabrier IMF, have emerged over 
the last decade to explain the observed distribution of stellar masses. The possibility of the universal nature of the stellar IMF 
pleads us to the interesting prospect of a universal mode of star-formation. It is well-known that turbulent fragmentation of gas in 
the interstellar medium produces a lognormal distribution of density which is further reflected by the mass-function for clumps at 
low and intermediate masses. Stars condense out of unstable clumps through a complex interplay between a number of dynamic 
processes which must be accounted for when tracing the origin of the stellar IMF. In the present work, applying the theory of 
gravitational fragmentation we first derive the mass function (MF) for clumps. Then a core mass function (CMF) is derived by 
> "\ allowing the clumps to fragment, having subjected each one to a random choice of gas temperature. Finally, the stellar IMF is 

O ! derived by applying a random core-to-star conversion efficiency, e, in the range of 5%-15% to each CMF. We obtain a power-law 

I \ 

O IMF that has exponents within the error-bars on the Kropua IMF. This derived IMF is preceded by a similar core mass function 

h ■ 

CZ2 which suggests, gravoturbulent fragmentation plays a key role in assembling necessary conditions that relate the two mass-functions. 
l—— 'In this sense the star-formation process, at least at low redshifts where gas cooling is efficient, is likely to be universal. We argue 

that the observed knee in the CMF and the stellar IMF may alternatively be interpreted in terms of the characteristic temperature 

^ ■ 

OO at which gas in potential star-forming clouds is likely to be found. Our results also show that turbulence in star-forming clouds is 

0\ ; 

'xT probably driven on large spatial scales with a power-spectrum steeper than Kolmogorov-type. 
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1. Introduction 

> : 

Explaining the origin of the stellar initial mass function 

b ■ 

C3 (IMF) is an outstanding question before the theory of star- 
formation. Ever since the seminal work by Edwin Salpeter 
(Salpeter 1955) on this subject, first suggesting a power-law 
mass-distribution of stellar masses for stars more massive than 
~ 1 M , astronomers have toiled hard to explain its origin. Early 
work on the subject attempted to explain it by invoking the pos- 
sibility of successive fragmentation of unstable clumps (Larson 
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1973). Though it was later shown that such fragmentation can- 
not continue unbridled and must be arrested at some point; the 
minimum mass of a fragment was calculated by Low & Lynden- 
Bell (1976). One of the first propositions for a universal IMF 
after Larson (1973), was made by Miller & Scalo (1979), who 
having studied a number of nearby star-forming regions pro- 
posed a lognormal distribution. The question has been receiv- 
ing greater attention in recent years, especially in view of the re- 
ported similarities between the distribution of core masses and 
the stellar IMF; see for example, Motte et al. 1998, Nutter & 
Ward-Thompson 2007, and Simpson et al. 2011, amongst a 
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number of other authors. 

While molecular clouds(MCs) contract under self-gravity, 
the presence of a turbulent velocity field fragments it into a 
number of clumps and filaments. With growing evidence point- 
ing towards the turbulent nature of MCs, a number of authors 
have attempted to simulate the effect of turbulence. Numer- 
ical models by for instance, Vazquez-Semadeni (1994), and 
Padoan & Nordlund (1999) showed that turbulence generates a 
lognormal density distribution in the volume of gas. Though 
turbulence is dissipative and decays relatively quickly, usu- 
ally on a timescale comparable to the sound-crossing time, it 
delays the onset of gravitational instability (see for e.g., Scalo 
& Pumphrey 1982, MacLow et al. 1998, Stone et al. 1998, and 
Klessen et al. 2000). Some of the more massive clumps once 
dominated by self-gravity, begin to collapse, where-as other 
smaller clumps are more likely to be sheared apart by turbulent 
flows. 

According to the fragmentation model introduced by Padoan 
& Nordlund (2002; PN), the mass distribution of clumps form- 
ing out of a turbulent density field reflects the distribution of 
gas density. In recent years two variants of this theory have 
been introduced. These theories have a subtle variation in the 
density threshold at which a volume of gas is likely to support 
fragmentation. In the PN theory this threshold is fixed by the 
thickness of the shocked layer resulting from collisions between 
turbulent flows. On the other hand, in the Krumholz & McKee 
(2005; KM) theory, the density threshold is set by equating the 
fastest growing unstable mode with the classical Jeans length. 
Finally, in the Hennebelle & Chabrier (2008; HQ theory, there 
is no density threshold but instead, density perturbations are as- 
sumed to condense out as unstable clumps. 

Numerical simulations by Padoan & Nordlund (2007) have 
demonstrated the power-law nature of clumps resulting from 
turbulent fragmentation. More recent simulations such as those 
by Federrath et al. (20 10), argue about the possible effects of 



the modes of driven turbulence, compressive against solenoidal. 
The resulting density distribution according to these latter au- 
thors depends on the driven mode. A possible short-coming 
of the theory of turbulent fragmentation is its failure to pre- 
dict beyond the mass-spectrum for gas clumps, and therefore 
the stellar IMF itself, since not all clumps actually spawn stars. 
While sophisticated numerical algorithms can possibly make 
the simulations more realistic by including effects of energy 
feedback, tracing the evolution of gas over several orders of 
magnitude, up to stellar densities, exacts enormous compu- 
tational costs. And so, a semi-analytic scheme like the one 
presented here can shed some light on the problem. 

Veltchev et al.(2Q\\) for example, recently presented one 
such calculation where the mass function for unstable clumps, 
cores, and stars was obtained. Having derived the mass function 
for clumps and cores using respectively, the turbulent fragment- 
ation theory, and the scaling relations due to Larson (1981), 
their work suffers on 2 counts : (i) the stellar masses were cal- 
culated by applying the Bondi-Hoyle(BH) accretion theory to 
collapsing cores. The BH theory as is well-known, assumes 
that the fluid being accreted originates from a non-self gravitat- 
ing background. This assumption breaks down in a collapsing 
core, (ii) The calculation of the velocity with which gas is 
accreted by the protostar is unclear. Despite this, their model 
provides a reasonable estimate of the protostellar accretion rate 
using which a power-law IMF has been derived. In a slightly 
earlier work (Anathpindika 2011), we had argued on similar 
lines where accounting for protostellar feedback, and treating 
the protostellar accretion rate as a free parameter, power- law 
IMFs were derived starting with different choices of core mass 
distribution that also included a uniform distribution. We had 
then sought to underline the importance of physical processes 
during the protostellar phase and the protostellar multiplicity in 
determining the final IMF. 

The problem of deriving the stellar IMF from first principles 



is indeed non-trivial and so a semi-analytic approach is use- 
ful to ascertain the effects of some crucial parameters. Cool- 
ing processes in star-forming clouds most likely hold the key 
to sub-fragmentation of unstable clumps into prestellar cores. 
However, owing to computational expenses, simulations about 
turbulent fragmentation of gas often assume a cold, isothermal 
gas, which though reasonable to some extent, fails to account 
for heating due to protostellar feedback and from other pos- 
sible external sources. Thus regions within a star-forming cloud 
exhibit substantial variation in gas temperature and therefore, 
an equally variable dynamical state. Possible effects of vari- 
ations in gas temperature were neglected by Veltchev et al. in 
their derivation of the core mass function. In the present work 
these effects are included by generating multiple realisations of 
clump fragmentation by subjecting them to a number of random 
choices of gas temperature. Having so derived a distribution of 
core masses, the stellar IMF is derived by applying a typical 
core-to-star conversion efficiency, e, between 5% - 15% to each 
bound core. Observations of star-forming clouds have yielded 
comparable values for the star-formation efficiency (e.g., Motte 
et al. 1998; Ward-Thompson et a/.2007). The paper is organ- 
ised as follows. The process of generating the mass-spectra is 
described in §2, and the spectra are then discussed in §3 and §4. 
We conclude in §5. 

2. Generating the mass functions 

We begin our calculations by assuming a volume of molecu- 
lar gas composed of the usual cosmic mixture where neutral 
Hydrogen makes up about 9 parts out of 10, while the remainder 
is made up by Helium and other elements like Carbon, Nitrogen 
and Oxygen. The average mass of a gas molecule, in, is taken to 
be 4x 1 0~ 24 gms. The gas within this volume is characterised by 
its average density, p avg ~ 10~ 21 gem -3 , temperature, T gas = 30 
K, and the magnitude of turbulent velocity, vo. The size, L, of 
the turbulent region was fixed using the velocity dispersion re- 



lation due to Larson (1981) so that the total mass within the 
volume, M tot ~ p aV gL 3 ~ 2.5 x 10 8 M . In a second realisation 
of the problem, the gas was maintained at a lower temperature, 
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15 K, wAM tot ~ 6xl0 5 Mq. The probability distribution 



of gas density, p(y), in the turbulent field is 
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where yo = -y, and 5 2 = ln(l + M 2 y 2 ) with y ~ 0.5, as sug- 
gested by Padoan & Nordlund (2002). M = v /a , is the Mach 
number of the turbulent field superposed on the gas maintained 
at a uniform temperature, T gas , and sound speed, oq. Finally, 
y = \n(p/p avg ). 

The turbulent field is defined by its energy-spectrum, &(k) oc 
k~ x , for x > 0, so that the turbulent velocity, v(k), is simply 



v 2 (k) = f &(k)dk = v\L x - x kl ~ X > 

Jk 



(2) 



where the integration and/or summation on the right-hand side 
extends over the driven wave-modes, k. Turbulence was driven 
over spatial scales between 10 pc and 200 pc, corresponding 
to wavenumbers between 1 and 8. The length of the fastest 
growing mode, Xf asU is given by the following equation, 

a r'ff\ u2 . m 
Afas '~(~Gp~) ' (3) 

p ~ M 2 p avg , where turbulence-induced shocks are assumed to 
be isothermal. The effective sound-speed, a e ff = («o + Vt), so 
that a slight manipulation of Eqn. (3) after substituting for p 
and a e ff, leads us to 



where Aj is the thermal Jeans length for sound-speed, ao. For 
a highly supersonic gas, the left-hand side of the above expres- 
sion may be approximated to 



A fast ~ —[kS(k)] 
v 
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or in terms of only the wave-vector as 
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The scale of fragmentation is set by Af ast , so that the mass of a 
clump, Mrfump, condensing out is of the order, 
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The number of clumps, N c i umps , condensing out is M, ,/M c i ump , 
so that the clump-mass function (MF), is 

dN c i um p S Mf f 



M 2 , 

clump 



k~ 3(i ~ x) . 



(8) 



dM dump 

The MF was generated by Monte-Carlo integrating the above 
equation with about 10 5 choices of random perturbations in the 
density field. We note that clump-formation was terminated 
once 2 Ndumps ^ M rot , so that mass is conserved. 

The dynamical stability of clumps, it is well-known, is de- 
termined by the balance between heating and cooling processes. 
Possible effects of various cooling processes such as molecu- 
lar line emission at densities below 10 4 crrT 3 , and the gas-dust 
coupling at much higher densities, were introduced in the cal- 
culations by subjecting the gas within these clumps to 10 3 ran- 
dom choices of temperature distributed uniformly between 7 K- 
20K, and then 6 K- 15 K, in each of the two realisations. The 
temperature distribution is characterised by a probability dens- 
ity, p(T), and the core mass function (CMF) was obtained by 
convolving it with the clump mass function, so that 



J p(T) dNd " mps ^ dT. 



(9) 



dM c i llmp: 

Alternatively, A core = f 22 ^ W as{ , so that M core ~ A 3 core p. The 
number of cores, N cor es, is then, M c [ ump /M core . Finally, the ini- 
tial mass function, (IMF), was derived by applying a random 
core-to-star conversion efficiency, e, between 5% to 15 %, to 
each bound core. 
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Figure 1 : Shown in each panel is the derived clump mass-spectrum for different 
choices of driven wave-modes, k, and the index of turbulent power-spectrum, 
x, as marked on each plot. The spectra are roughly lognormal in nature, though 
with a slight power-law extension at the high-mass end. 



3. Mass functions (MF) 



Using the basic turbulent fragmentation hypothesis we have 
derived mass functions (MFs) for dense clumps, cores and stars. 

4 



Fragmentation of the volume of gas considered here is purely 
due to the growth of density perturbations that are random in 
nature. While our interest lies in the clump-formation via frag- 
mentation, the dissipation of turbulent energy is neglected. Res- 
ults of this exercise are summarised by the plots in Figs. 1 
through to 5. The distribution of clump masses derived here has 
been shown in Fig. 1 . Irrespective of the choice of driven wave- 
modes, the mass-functions exhibit remarkable similarities. The 
distributions in each case are roughly lognormal, though with a 
slight power-law extension at the high-mass end. The spectra, 
however, are significantly steeper than that for typical clumps, 
which have a slope of ~ 2 (e.g. Kramer et al. 1998). Also, 
it follows from Eqn. (5), a relatively shallow energy spectrum 
(x=1.2) shows a greater tendency of producing larger clumps. 
This can indeed be seen from the plots in Fig. 1 

The core mass function(CMF), generated for different 
choices of the spectral index, x, and driven modes, k, have been 
shown in Fig. 2. Overlaid on the plot shown in the upper panel 
of Fig. 2 is the power-law CMF defined as, 
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< 1, 



0.02 <fr< 0.1. 



The power-law with these exponents fits the derived CMF reas- 
onably well, acceptable at 5% significance level of the^f 2 -test, 
performed with 5 degrees of freedom, viz., the gas temperature, 
density perturbations, turbulence Mach number, M, spectral- 
index, x, and the wavemodes driven, k. 

Evidently the power-law defined by Eqn. (10) is consider- 
ably steeper than that for cores found in typical star-forming 
regions such as the Orion MC (e.g., Nutter & Ward-Thompson 
2007), and the Ophiuchus MC (e.g., Johnstone et al. 2001). In- 
creasing the gas temperature and the Mach number, Al, raises 
the fragmentation length and therefore the mass of resulting 



cores which results in a marginally shallow CMF (M^ r ^ s , for 
M cores /M e > 1), as can be seen in the central panel of Fig. 
2. Likewise, a shallow CMF also results from driving longer 
wavemodes, k = (6,8), when comparatively lesser power 
reaches smaller spatial scales. This latter CMF is shown in the 
lower panel of Fig. 2. A unique feature of the CMFs derived 
here is the appearance of a plateau for intermediate masses in 
the range (0.2,1) M , which is probably due to the uniform dis- 
tribution of temperature used in generating the CMF. 

The core-to-star conversion efficiency, e, is the sixth para- 
meter in this scheme apart from the five listed previously. The 
derived IMF has been shown in the two panels of Fig. 3, where 
the one in the upper panel was derived from the CMF defined by 
Eqn. (10). It is significantly steep for M star /M Q > 0.5. On the 
other hand, the IMF derived from the relatively shallow CMF 
obtained by driving the k = (6, 8) modes, has been shown in the 
lower-panel of Fig. 3. It has been fitted by a power-law defined 
below, 
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0-01 < < 0.04, 



acceptable at a 10% significance level of the ;^ 2 -test with 
6-degrees of freedom. More interestingly, the power-law 
exponents are within the error-bars on the universal IMF 
suggested by Kroupa (2002). It can be seen that like the CMF, 
this IMF also has a plateau in the sub-Solar mass range which 
most likely reflects a similar feature of the CMF. 

Lowering the initial gas temperature (T gas = 15 K, 
M=15) 

Retaining all other parameters, the mass function for clumps 
and cores was re-derived in this case, the result of which has 
been plotted in Figs. 4 and 5. Lowering the initial temperature 
of the gas, T gas , by a factor of 2 to 15 K, does not much alter the 
overall nature of the two distributions, though massive clumps 








x=1.2 






x=3.67 




1 


k=(6,8);M=l5 

\ 


; 1 













Figure 2: Shown in the three panels are the core mass functions derived for 
different choices of driven modes, k, and the power- spectrum index, x. Plotted 
with a continuous black line are the power-laws discussed in the text above. 



Figure 3: The stellar initial mass function (IMF) derived for two choices of 
k has been shown in the two panels here. The power-law fit to the derived 
spectrum in the lower panel is defined by Eqn. (1 1) in the main text. 

with mass in excess of 10 4 M now do not form as can be seen 
in Fig. 4. As before, gas within these clumps was subjected 
to a random choice of temperatures in the range 6 K- 15K, 
and then allowed to fragment. The resulting distribution of 
core masses has been shown in the upper-panel of Fig. 5. The 
distribution of core-masses in the lower-panel was derived for a 
slightly higher initial Mach number (M = 20). Not only is the 
CMF in this latter case relatively shallow (oc M^ s ; M cores > 1 
M ), but the plateau at intermediate masses is also less wide. 

4. Discussion 



In the previous section, assuming a set of demonstrative 
parameters we derived mass distributions for clumps, prestellar 
cores and stars by invoking the theory of turbulent fragment- 



x-1.2 - 
x-1.67 
x=3.67 



k=(6,8), Mach=15 




Figure 4: The clump-mass spectrum derived in the next case where the gas 
was originally maintained at T gas = 15 K. Apart from a broader full-width at 
half-maximum, these spectra do not differ much from those plotted in the upper 
panel. 

ation. Analogous to the PN02 model, we commenced with a 
lognormal PDF while setting the density threshold as dictated 
by the jump condition for an isothermal shock. The derived 
clump mass function largely reflects the density PDF, though 
it also has a power-law tail at the high-mass end. This tail 
is relatively shallow in the case with T gas = 15 K. Steepening 
the power-spectrum index, on the other hand, simply shifts the 
mass distributions towards lower masses, closer to those of typ- 
ical clumps. 

The core mass function, irrespective of the choice of 
wavenumbers driven, the temperature of gas within clumps or 
indeed, the Mach number, has two parts : a Salpeter-like power- 
law for M core /M Q > 1, and a plateau at lower masses. We have 
also seen that the power-law for M core /M > 1 is relatively 
steep, and lies between 3.5-4 which is probably the result of 
an analogously steep mass-function for clumps which precede 
the CMFs. Interestingly though, these CMFs extend into the 
regime of low mass cores (< 10~' M ). The plateau observed 
in the CMFs also appears in the stellar IMFs as can be eas- 
ily noticed through a comparison between Figs. 2 and 3. The 
CMFs and the IMFs derived here both show a significant over- 
population in the (0.5,1.0) M range. The plateau in the CMFs 
derived here merely reflect the distribution, p(T), of gas tem- 



perature within clumps which, we remind, was assumed to be 
distributed uniformly between 7 K - 30 K, and then between 6 
K- 15 K in the next realisation. The plateau can be avoided by 
confining the gas temperature to a fairly narrow range implying 
that, the temperature of gas across star-forming clumps is un- 
likely to show much variation. Of course, we are discounting 
regions which may exhibit extreme conditions. 

Fragmentation in a volume of turbulent gas proceeds along 
a hierarchy that begins with clumps and filamentary clouds. 
Some of which fragment further to produce prestellar cores. 
Indeed, some recent maps of star-forming regions in various 
sub-millimetre wavebands have revealed extensive networks of 
filamentary clouds with a number of prestellar cores (e.g. Andre 
et al. 2010; Men'shchikov et al. 2010). From the preced- 
ing analysis in §2 it is clear that the prevalent physical condi- 
tions within gas clumps are crucial in determining the procliv- 
ity towards further fragmentation. To the first order, the min- 
imum mass of fragmentation is set by the thermal Jeans mass, 
My ~ pA 3 ,. In reality of course, the magnetic field will revise 
this mass upward, however, in the present exercise we do not 
include the effects of magnetic field. 

Gas temperature and density are locked in a reciprocally 
causative relationship. It is in this respect that a density 
threshold assumes critical significance, for cooling mechanisms 
are density dependent. Atomic cooling for instance, dominates 
in the relatively low-density gas up-to a few hundred particles 
per cubic centimetre, then the coupling between gas and dust 
dominates up-to densities of a few thousand particles per cubic 
centimetre and at even higher densities such as those in typical 
cores (> 10 5 cirT 3 ), molecular species freeze out. Similarly, ef- 
fects of magnetic flux diffusion are likely to become significant 
at densities upward of several thousand particles per cubic cen- 
timetre. Taken cognitively, these effects can raise the minimum 
mass of fragmentation considerably (e.g. Larson 2005). 

Here we have only included the effects of gas thermodynam- 
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Figure 5: Shown in the upper and lower-panel of this figure are the CMFs 
derived in the second realisation of the problem where the original gas was 
maintained at 15 K. A power-law fitted to the derived distribution has been 
plotted on top with a black line. 



ics by accounting for the gas temperature and the velocity dis- 
persion. The contribution of the velocity dispersion to the frag- 
mentation length scale is determined by the energy-spectrum, 
S(k). In general, the impact of the energy-spectrum is readily 
visible on the mass-function for clumps where relatively smal- 
ler clumps form with steepening energy-spectrum; see plots in 
Fig. 1. The CMFs on the other hand, while exhibiting im- 
munity to changes in S(k), are significantly affected by changes 
in the thermal state of the gas, as demonstrated by plots shown 
in the central panel of Fig. 2. This is further corroborated 
by results from the exercise where the temperature of the pre- 
fragmentation gas was lowered by a factor of 2; see also Figs. 4 
and 5. In either case, temperature of the gas within a clump was 
treated as a random variable to study the effects of its variation. 
The gas temperature thus, was allowed to take values within 
a range typically found in star-forming clouds such as Taurus 
(e.g., Falgarone et al. 1998), Ophiuchus (e.g., Johnstone et al. 
2000), and Orion (e.g., Yu et al. 1997). 

Prestellar cores collapse under favourable conditions to form 
stars and dwarfs, though the collapse phase itself is associated 
with competing processes such as feedback from young stellar 
objects, competition between accreting stars to acquire mass 
and the interaction between a young star and its attendant disk. 
These processes, though highly non-linear in nature, are crit- 
ical in determining the stellar IMF. In literature, these processes 
have sometimes been grouped together as the core fragmenta- 
tion mass function (e.g. Elmegreen 2000; 201 1). Clearly, stellar 
multiplicity is an essential component of the stellar IMF (e.g., 
Goodwin & Kouwenhowen 2009). In an earlier work (Anath- 
pindika 201 1), we had presented a model of core fragmentation 
that included the details listed above, and therefore the prospect 
of forming brown-dwarfs via core-fragmentation. In view of the 
non-triviality about estimating the mass of a protostellar disk, 
that model failed to predict any objects possibly resulting from 
disk-fragmentation. However, assuming the rate of protostellar 



accretion as a free parameter, we had then demonstrated that the 
processes during the collapse of a core hold the key to determ- 
ining the nature of the IMF. That model had also demonstrated 
the robustness of the IMF to changes in the preceding CMF 
In the present work though, we have derived the IMF by apply- 
ing a typical core-to-star conversion efficiency, e, to each bound 
core. In the present case it is therefore not possible to comment 
on the companion stars. However, it is clear that the power- 
law tail in the CMF reappears in the IMF, albeit with a slightly 
shallower slope and with significant changes to the plateau at 
sub-Solar masses. 

The number of parameters used in the current exercise could 
invite some criticism, but these parameters are essential com- 
ponents of the turbulent fragmentation model and adopted from 
recent reports about star-forming clouds. For a typical choice of 
parameter values, the derived mass functions are broadly sim- 
ilar to those determined observationally. The agreement was 
best seen for results obtained with - (i) a steep power-spectrum, 
typically steeper than the Kolmogorov-type, up to about k~ 3 61 , 
and (ii) by driving longer wavemodes, k = (6, 8). In view of 
(i), we suggest that the power-spectrum index for gas in MCs 
could indeed be steeper than Kolmogorov-type. Second, turbu- 
lence is likely to be driven on larger scales which is consistent 
with observations and has also demonstrated by Brunt et al. 
(2009) through hydrodynamic simulations. Furthermore, stud- 
ies of the 13 CO emission line for the Cygnus star-forming re- 
gion (Schneider et al. 201 1), and a number of other star-forming 
clouds in the Galactic ring survey (Roman-Duval et al. 201 1), 
have shown that turbulence is indeed driven on larger spatial 
scales, typically on the scale of the the cloud size. A point that 
has not been addressed in the present work is about the likely 
effects of changes in the gas density. However, from Eqns. (3) 
- (6) it is quite clear that lower density may very well produce 
predominantly larger clumps, which might possibly explain the 
relatively steep CMFs observed in low-density clouds. It might 



also explain the steep stellar mass distributions in some spiral 
galaxies with peculiar properties, such as those with low surface 
brightness (e.g. Lee et al. 2004). Thus, if indeed the stellar IMF 
has a universal nature then the underlying processes which gen- 
erate it could also be similar. Gravoturbulent fragmentation, it 
appears, could play a key role in assembling the initial condi- 
tions for star-formation. This theory also admits the possibility 
of rapid star-formation which makes it more attractive. 

5. Conclusions 

Using the concept of turbulent fragmentation we have de- 
rived the mass spectra for clumps, prestellar cores and stars via 
a Monte Carlo integration of the underlying probability distri- 
bution. We first calculated the clump mass spectrum by al- 
lowing the random density perturbations to condense out as 
clumps, while the CMF was derived through the convolution of 
this mass spectrum with a random distribution of gas temper- 
ature for the dense clumps. Irrespective of the initial gas tem- 
perature, we observed that best agreement with observational 
findings resulted from driving longer wavemodes and with a 
power-spectrum steeper than Kolmogorov-type. Thus, while 
fragmentation of the original volume of gas produces clumps 
which themselves may break-up further, chemical processes 
that determine gas temperature in the dense phase are likely to 
play a crucial part, leading up to the distribution of core masses. 
In fact, we also suggest that the characteristic knee seen in the 
CMF and the stellar IMF could be a pointer to the temperature 
of the gas within star-forming clumps, and therefore, the effi- 
ciency of cooling processes. For the theory of fragmentation to 
predict the core masses more accurately, particularly at the low- 
mass end, it must be coupled with appropriate thermodynamic 
details of the density field. Arguably, this is a non-trivial task, 
and so, in the present work temperature has been substituted as 
a random variable that takes values similar to those found in 
typical star-forming clouds. 



The IMF was also derived in a similar fashion by treating 
the star-formation efficiency as a random variable. The IMF 
so derived was found to be in reasonably good agreement with 
the universal IMF, acceptable at a 10% significance level of the 
^ 2 -test with 6 degrees of freedom, viz. the gas temperature, 
density perturbations, the Mach number and spectral-index of 
the turbulent field, the wavemodes driven, and the core-to-star 
conversion efficiency. It is therefore clear that apart from the 
physical properties of the gas, the history of protostellar col- 
lapse is critical towards determining the nature of the IMF. Our 
work shows that a typical CMF, and therefore an IMF, emerge 
quite naturally from the application of the theory of turbulent 
fragmentation by driving longer wavemodes. This claim is cor- 
roborated by numerical and observational results presented by 
Brunt etal. (2009) and Schneider et al. (2011). While isol- 
ated regions within star-forming clouds may exhibit extreme 
physical conditions leading to corresponding variations in the 
IMF, results derived here suggest that most star-formation in 
the nearby clouds is unlikely to vary much; see also the review 
by Bastian e/a/.(2010). 
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